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ABSTRACT 

X-ray satellites since Einstein have empirically established that the X-ray luminos- 
ity from single 0-stars scales linearly with bolometric ImTiinosity, Lx ~ 10~''iboi- But 
straightforward forms of the most favored model, in which X-rays arise from instability- 
generated shocks embedded in the stellar wind, predict a steeper scaling, either with 
mass loss rate ^ M ^ Lj^Jj if the shocks are radiative, or with ~ ~ i^Qi 
if they are adiabatic. This paper presents a generalized formalism that bridges these 
radiative vs. adiabatic limits in terms of the ratio of the shock cooling length to the 
local radius. Noting that the thin-shell instability of radiative shocks should lead to 
extensive mixing of hot and cool material, we propose that the associated softening 
and weakening of the X-ray emission can be parameterized as scaling with the cooling 
length ratio raised to a power m, the "mixing exponent." For physically reasonable 
values TO « 0.4, this leads to an X-ray luminosity ^ M°'^ ^ Lboi that matches 
the empirical scaling. To fit observed X-ray line profiles, we find such radiative-shock- 
mixing models require the number of shocks to drop sharply above the initial shock 
onset radius. This in turn implies that the X-ray luminosity should saturate and even 
decrease for optically thick winds with very high mass-loss rates. In the opposite limit 
of adiabatic shocks in low-density winds (e.g., from B-stars), the X-ray luminosity 
should drop steeply with Af^. Future numerical simulation studies will be needed to 
test the general thin-shell mixing ansatz for X-ray emission. 

Key words: shock waves - stars: early-type - stars: winds - stars: mass loss - X-rays: 
stars 



1 INTRODUCTION 

Since the 1970's X-ray satellite missions like Emstem, 
ROSAT, and most recently Chandra and XMM-Newton 
have found hot, luminous, O-type stars to be sources of 
soft (< 1 keV) X-rays, with a roughljQ linear scaling 
between the X-ray luminosity a nd the stellar bo l omet- 
ric luminosity, ix_~ lO'^Lboi jChlebowski et al.l 19891: 
Kudritzki et all 119 96*: 'Bcrghocfcr c t all 1 19971: ISana et'al] 
20061 : iGiidel fc Naza .2009 : .Naze .20091 : iNaze et all l201lD . 
In some systems with harder (a few keV) spectra and/or 
higher L^, the observed X-rays have been associated with 
shock emission in c o lliding wind binary (CWB) systems 
IStevens et al.l Il992l : iGagnel l201l|). or with magnetically 
confined wind shocks fMCWS) (|Babel fc Montmerld[l997l : 



i.e., extending over 2 dex in Ltoi, v^ith a typical scatter of 
' ±0.5 dex 



IWaddl201ll ). But in putatively single, non-magnetic O-stars, 
the most favored model is that the X-rays are emitted from 
embedded wind shocks that form from the strong, intrinsic 
instability (the "line-deshadowing instability" or LDI) asso- 
ciated with the driving of the se winds by line-scatte ring of 
the s t ar's radiative fiux ( Owock i ct al. 1988; Fcldnici cr et al.l 



ll997l : lDessart fc Owockill2003l : ISundgvist fc Owockill2012l ). 



This LDI can be simply viewed as causing some small 
(< 10~^) fraction of the wind material to pass through an X- 
ray emitting shock, implying in the case that the full shock 
energy is suddenly radiated away that the X-ray luminos- 
ity should scale with t he wind mas s loss rate, Lx ~ M. 
But within the standard lCastor et al. 1 (|1975| . hereafter CAK) 
model for such radiatively driven stellar winds, this mass loss 
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Figure 1. Summary sketch of the key concepts and results of the scaUng analysis in this paper. The illustration of the cooling zone 
from a wind shock shows associated scalings for X-ray luminosity Lx with mass loss rate M and bolometric luminosity Lbob based on 
the conversion factor of wind kinetic energy into X-rays, which for radiative shocks is a constant, but for adiabatic shocks is reduced 
by the ratio of the radius to cooling length, r/£ <^ 1. (See eqn. 1121 in text.) In addition, thin-shell mixing of such radiative shocks is 
then posited to lead to a reduction of the X-ray emitting fraction that scales as a power-law of the cooling length, fcx ~ For CAK 
wind index a, a mixing exponent m = I — a leads to the observationally inferred linear scaling of X-rays with bolometric luminosity. 



rate increases with luminositjjfl as M ~ L^/j^ 
the latter scaling uses a typical CAK power index a w 0.6 
(jPuls et al.l 120001 1. This then implies a super-linear scaling 
for X-ray to bolometric luminosity, Lx ~ ^boij which is too 
steep to match the observed, roughly linear law. 

In fact, the above sudden-emission scaling effectively as- 
sumes the shocks are radiative, with a cooling length that 
is much smaller than the local radius, £ ^ r. In the op- 
posite limit i ^ r, applicable to lower-density winds for 
which shocks cool by adiabatic expansion, the shock emis- 
sion scales with the X-ray emission measure, EM~ J p^dV, 
leading then to an even steeper scaling of X-ray vs. bolo- 
metric luminosity, Lx ~ M'^ ~ ^boi- 

Both these scalings ignore the effect of bound-free ab- 
sorption of X-rays by the cool, uns hocked material that rep- 
resents the bulk of the stellar wind. lOwocki fc CohenI (|l999l . 
hereafter OC99) showed that accounting for wind absorption 
can lead to an observed X-ray luminosity that scales linearly 
with Lboi, but this requires specifying ad hoc a fixed radial 
decline in the volume filling factor for X-ray emitting gas. 
We show below ( H2.4|l that this filling factor should actually 
be strongly affected by the level of radiative cooling. More- 
over, while modern observations of spectrally resolved X-ray 
emission profiles by Chandra and XMM-Newton do indeed 
show the expected broadening from shocks embedded in the 
expanding stellar wind, the relatively weak blue-red asym- 
metry indicates that absorption effects are modest in even 



^ For simplicity, this ignores a secondary scaling of luminosity 
with mass; see i]3.4l 



the densest OB-star winds (|Cohen et al.ll2010l . [ioTH ). Since 
many stars following the I/x-Lboi empirical law have weaker 
wind s that are largely optically thin to X-rays (jCohen et al.l 
ll997l : lNaze et aLlbOllh . it now seems clear that absorption 
cannot explain this broad Lx scaling. 

The analysis here examines instead the role of 
radiative cooling, and associated th i n-shell instabilities 
(IVishniadll994l: IWalder fc Folinil Il998l : ISchure et al.] l2009l : 
iParkin fc PittardI bOlof T in mixing shock-heated material 
with cooler gas, and thereby reducing and softening the over- 
all X-ray emission. As summarized in figure 1, for a simple 
parameterization that this mixing reduction scales with a 
power (the "mixing exponent" m) of the cooling length, , 
we find that the linear Lx-Lboi law can be reproduced by 
assuming m ~ 0.4. The development below quantifies and 
extends a prelimina ry conference presen tation of this thin- 
shell mixing ansatz jOwocki et aLlfioill ). 

Specifically, to provide a quantitative basis for bridging 
the transition between radiative and adiabatic shock cool- 
ing, the next section (§2) first analyzes in detail the X-ray 
emission from a simplified model of a single, standing shock 
in steady, spherically expanding outflow. The following sec- 
tion (§3) then generalizes the resulting simple bridging law 
to account for thin-shell mixing, and applies this in a sim- 
ple exospheric model to derive general scalings for X-ray 
luminosity from a wind with multiple, instability-generated 
shocks assumed to have a power-law number distribution 
in wind radius. A further application to computation of X- 
ray line profiles (§4) provides constraints on the mixing and 
shock-number exponents needed to match observed X-ray 
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emission lines. Following a brief presentation (§5) of full in- 
tegral solutions for X-ray luminosity to complement the gen- 
eral scaling laws in §3, the final section (§6) concludes with 
a brief summary and outlook for future work. 



2 STANDING SHOCK MODEL 
2.1 Energy Balance 

Most previous analy s es of X-rays from massiv e stars (e.g., 
IWoidowski fc Schulj l2005l : ICohen et all I2OIII ) have been 
cast in terms of a density-squared emission measure from 
some fixed v olume filling factor for X-ray emitting gas 
(see however ^Krolik fc Raymond 19851: iHillier et all 1 19931 : 
iFeldmeier eTal . 1997; Antokhin et all 120041 '). So let us be- 
gin by showing explicitly how that picture must be modified 
to account for the effects of radiative cooling, which can be 
important and even dominant for O-star wind shockfl By 
focusing on the simple example of a steady, standing wind- 
shock, it is possible to carry out an analytic analysis that 
derives, more or less from first principles, a simple bridging 
law between the scalings for adiabatic vs. radiative shocks. 

Specifically, let us consider an idealized model in which 
a spherically symmetric, steady-state stellar wind with mass 
loss rate M and constant, highly supersonic flow speed 
Voo undergoes a strong, standing shock at some fixed ra- 
dius r — Vs. Relative to an unshocked wind with density 
Pw ~ M/AnVoor^, the post-shock flow at r > rs has a den- 
sity that is a factor p/pw ~ 4 higher, with a post-shock 
speed that is a factor v/Voo = 1/4 lower. The reduction in 
flow kinetic energy results in a high, immediate post-shock 
temperature. 



3 I^V^ 
16 k 



= 14 MK 



1000 kms' 



(1) 



where k is Boltzmann's constant, and the latter evaluation 
assumes a standard molecular weight /j, — 0.62mp, with TUp 
the proton mass. 

But following this sudden shock increase, the combined 
effects of adiabatic expansion and radiative cooling cause 
the flow temperature T to decrease outward. For pressure 
P — pkT/fj. and internal energy density e = (3/2)P, the 
steady-state energy balance for a general vector velocity v 
is 



V • (ev) 



-PV- v-p^An,(r), 



(2) 



where Am = A/(^,ofi^), with A(r ) the optically thin cooling 
function fe.g.. lSmith et aLll200ll ). and pe and pp respectively 
the mean mass per electron and per proton. 

In general, we should also include a detailed momentum 
equation to account for possible acceleration of the post- 
shock flow, for example from the inward pull of stellar grav- 
ity, or the outward push of the gas pressure gradient. But 
the analysis here is greatly simplified if we make the reason- 
able assumption that these two countervailing accelerations 
roughly cancel, and so give a constant speed v = Voo/4 for all 



^ This was actually noted explicitly by IZhekov fc Fallal |20o3), 
but their results were nonetheless still cast in terms of a density- 
squared emission measure that is not appropriate for radiative 
shocks. 



r > Ts. For this case of a steady, spherical, constant-speed 
post-shock outflow, the vector energy equation ((2]) reduces 
to a simple differential equation for the decline in tempera- 
ture with radius r. 



dT 
dr 
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'3 7 



2p 

3kv 



pAm(T) , 



(3) 



wherein the first and second terms on the right-hand-side 
respectively represent the effects of adiabatic expansion and 
radiative cooling. This can alternatively be cast in terms of 
a temperature scale length. 



1 _ 1 dT _ 4 Kcp 



where 



Atc(r) 



8pAn,(r) 

kTVoc 



(4) 



(5) 



is a mass cooling coefficient (with CGS units cm^g"^), de- 
fined as the inverse of a characteristic cooling column mass. 
The corresponding cooling length is given by i = 4/kcP = 
1/^cPw, defined such that the radiative and adiabatic cool- 
ing terms are equaQ when I = r. 



2.2 X-ray luminosity 

For any local post-shock temperature T, let f-^{T) represent 
the fraction of radiation emitted in an X-ray bandpass of 
interest. Neglecting for now any wind absorption, the to- 
tal X-ray luminosity from this single standing shock is then 
given by radial integration of the associated X-ray emission, 



47r / p'A„,(r)/x(T)r'dr. 



(6) 

Using Q, this can be recast as an integral over temperature, 
' 3 p'Am(T) 



Lxs = 127r 
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(7a) 
(7b) 



where the latter approximatior(f| uses single-point trape- 
zoidal integration, assuming that /x declines from its post- 
shock value fx{Ts) to zero over a temperature range 5Ts 
from the initial post-shock temperature Ts = T{rs). Here 
Pws = M/{A-KVoors) is the wind density just beforethe shock. 



* The ratio £/r differs only by an order-unity factor from the ratio 
of coo ling to escape time, x = tcool/*esc, defined by I Stevens et all 
to characterize the transition from radiative to adiabatic 
shocks in colliding stellar winds. 

^ Aside from /x(T), the combination of other factors in the in- 
tegrand for I l7al l becomes constant in T in the radiative limit 
Kcpr 2> 1, and scales as T^"^/** in the opposite, adiabatic limit. 
For /x(r) that declines roughly linearly with temperature, the ap- 
proximate trapezoidal integration l|7b|l thus becomes nearly ex- 
act in the radiative limit, while mildly underestimating the ac- 
tual value in the adiabatic limit, for example by about 15% for 
5Ts = Ts/2. 
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(8a) 
(8b) 
(8c) 



with Tj = Ts/l(fK and Tkev = kT^/keV. The numerical 
evaluation in H8c|) assumes an approximate fit to the cooling 
function, A(rs) ^ 4.4 x 10"^^/^^ ergcm^s"^ over the rel- 
evant range of shoc k temperatures, 10^'^ K < Tg < K 
ijSchure et al.ll2009l ^. Recall that the shock cooling length is 
set by £b = 1/kcsPws- 

For context, the mass absorption coefficient for bound- 
free absorption of X-rays when smoothed over bound-free 
edges, also roughly follows an inverse-square scaling with 
energy. Over the relevant energy ran ge 0.5 — 2 keV, w e 
can use the the opacity curves of. e.g..iCohen et all (|2010t ). 
iLeutenegger et alj (|2010t ). or lHerve et al.l (|2012l ) to write an 
approximate scaling form. 



30- 



(9) 



where i?kcv is now the X-ray photon energy in keV. 

By casting the cooling strength in an opacity form nor- 
mally used to describe absorption, we are thus able to make 
direct comparisons between cooling and absorption, and so 
characterize their respective domains of importance. This is 
further facilitated by the quite fortunate coincidence that 
both have similar inverse-square scalings with their associ- 
ated energy. 

Since these respective energies are usually roughly com- 
parable, -Ekev ~ 2kev, the fact that the numerical factor for 
Kcs is about 25 times greater than for Kbf means that cool- 
ing can become important even in winds that are too low 
density to have significant absorption. As such, winds with 
adiabatic shocks are always optically thin, whereas shocks 
in optically thick winds are always radiative. Moreover, as 
detailed in even in the bulk of O-star winds for which 
X-ray absorption is weak or marginal, the structure of X- 
ray emission associated wind shocks should be dominated 
by radiative cooling. 



2.3 Bridging law between radiative and adiabatic 
limits 

Noting that the pre-shock kinetic energy luminosity of the 
wind is 



L^ = 2-Kr!p„,V^ = MV^/2, 



(10) 



we can use (|8ap to eliminate the cooling function A,n from 
l|7bp , recasting the X-ray luminosity scaling as a "bridging 
law" between the radiative and adiabatic shock limits, 



16 1 + ls/r, 



(11) 



where /xs = fx{Ts)5Ts/2Ts is now a cooling-integrated shock 
X-ray fraction. 

Note here that the combination of factors multiplying 



Lu) on the right-hand-side of (|lip is just the kinetic energy 
conversion factor introduced in the summary figure [Jl 



kv — — 



/x; 



16 1 + £s/r, ■ 
For high-density, radiative shocks with is <^ rs 
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16 



/x 



(12) 



(13) 



As a physical interpretation, 9/16 is just the fraction of wind 
kinetic energy that is converted to post-shock heat, which is 
radiated away before any losses to adiabatic expansion, with 
the fraction /xs emitted in the X-ray bandpass of interest. 
Note that this scales linearly with density and thus mass loss 
rate, showing that a standard density-squared emission mea- 
sure does not represent an appropriate scaling for emission 
from radiative shocks. 

For lower- density, adiabatic shocks with 3> ^s, this X- 
ray emission is reduced by the ratio Vs/is, giving the scalings. 



/x: 



_9_ 

16 



— /xs g ' s '^cxy^csPws 

= /xs 487rrfA,^(rs)pL 
3An.(rs) ^M 



/x; 



(14a) 

(14b) 
(14c) 

(14d) 



which thus recovers the density-squared scaling for X-ray 
luminosity, showing that emission measure does provide an 
appropriate scaling for emission from adiabatic shocks. 

Finally, note that the general bridging law (|ll|l can al- 
ternatively be written as a modification of either the radia- 
tive or adiabatic scaling, 

-^xm,rad 



Lx 



1 + 4/rs 

-^'xm,ad 



(15a) 
(15b) 



1 + rs/4 

This shows that the X-ray luminosity is always limited to 
be somewhat below the smaller of the radiative or adiabatic 
luminosities, i.e. Lxs ^ min(Lxm,rad, ixm,ad). 

2.4 Local X-ray emissivity 

To facilitate application of these single standing-shock scal- 
ings to the more complex case of multiple embedded wind 
shocks generated from the line-deshadowing instability, let 
us next recast these results in terms of the local X-ray emis- 
sion from an individual shock. Since X-ray emission arises 
from collision of ions and electrons, it is common to write 
the X-ray emissivity (per unit volume and solid angle) as 
scaling with the square of the local density. 



»7x = CsfvP 



(16) 



where Cs is a constant that depends on the shock strength 
and atomic physics, and /v represents a local volume filling 
factor for shocked gas that is sufficiently hot to emit X-rays. 
If each individual shock has an associated filling factor /vs 
proportional to its post-shock cooling length, then the total 
filling factor from an ensemble of shocks can be written, 

dNs 



fvir) = /vs 



dlnr 



= /vs ns(r-) , 



(17) 



Thin-shell mixing and 0-star X-rays 5 



where Ns{r) is the cumulative number of shocks up to ra- 
dius r, and Us measures the local differential number of 
new, emerging shocks. This formalism emphasizes the im- 
portance of the number of shocks and their spatial distri- 
bution; compared to the traditional emission measure ap- 
proach, it should provide more physically motivated con- 
straints for wind-shock X-ray production in massive stars. 

In particular, the analysis below of multiple, instability- 
generated shocks assumes a power-law scaling for rig (see eq. 
[25]); but for the above single-shock model, this just takes 
the form of a Dirac delta-function, Us = rS{r — rg). The 
associated X-ray luminosity is then given by integration of 
the emissivity rj^ over solid angle and volume. 



16n 



Cs/vs rS{r — rs)p^ r^dr . 



(18) 



Upon trivial evaluation over the delta function, comparison 
with the scalings (|14c|) and ()15b|) yields the identifications, 



and 



Cs — — Am(T's)/xi 



/v: 



1 + 



(19) 



(20) 



This thus now sets a bridging law at the level of an indi- 
vidual shock, with /vs characterizing the fraction of single- 
shock emission measure that actually contributes to radia- 
tive emission. The adiabatic limit 1^ ^ t's gives /vs ~ 1 and 
so a density-squared scaling for the emissivity (|16|l : the ra- 
diative limit reduces this by the small factor Is/ts <^ 1, thus 
giving this emissivity a linear-density scaling set by the total 
kinetic energy flux through the shock. 

Note that the factor Cs depends only on the shock 
strength (through the post-shock temperature Tg), while the 
shock filling factor /vs depends on the local shock radius and 
cooling length, as set by the pre-shock wind density. Thus, 
although these scalings are derived from a simple model of a 
single, standing shock, they should also be generally appli- 
cable to more complex, moving shock structures, if specified 
in terms of the pre-shock wind densitjjf] and the relative ve- 
locity jump that sets the shock strength and temperature. 



3 X-RAYS FROM INSTABILITY-GENERATED 
EMBEDDED WIND SHOCKS 

3.1 Exospheric scaling for Lx 

To model the X-ray emission from multiple, embedded wind 
shocks, let us next write the X-ray luminosity in a fully 
general, 3-D form that accounts for possible directional de- 
pendencies in emission and absorption, 

ix = y"d^-y"dr!77x(r,n) (21) 



^ For instability-generated wind structure. X-rays can arise from 
collisions between clumps that have been comp ressed to some 
fraction of the wind volume l lFeldmeier et al .Il997l) . implying then 
a higher input density that would lower the cooling length. For 
simplicity, the analysis here does not account for this possibility, 
but it could enhance the importance of radiative cooling, leading 
to an even larger effective ratio Kcs/^bf of cooling to absorption. 



where the optical depth t{v, n) accounts for bound-free ab- 
sorption of X-rays emitted at location r in the observer di- 
rection n. 

In principle, instability-generated. X-ray emitting 
shocks will be associated with a comp lex, 3D, stochastic 
wind structure (|Dessart fc Owockill2003l l . But upon averag- 
ing over small scales, this can be reasonably well described 
by a globally spherical wind emission model (OC99). In 
modeling line emission, accounting for the observed Doppler 
shift from wind expansion still requires including a direc- 
tional dependence of the emis sivity and optical depth; see 
^and lOwocki fc CohenI (|200ll . hereafter OCOl). But to de- 
rive the general scalings for the X-ray luminosity, one can 
again take the total X-ray emissivity to be an isotropic func- 
tion of the local radius, ri^{r). Moreover, given the weak to 
moderate importance of absorption in all but the densest 
winds, its overall role in scaling relations can be roughly 
taken into a ccount through a simple exospheric approxima- 
tion (OC99. rLeutenegger et aLlbOlol l. for which the integra- 
tions over solid angle and volume in ^21\ reduce to just a 
single integration in radius. 



Lx 



167r' 



'7x('") dr 



(22) 



where the integral lower bound is taken from the larger of 
the X-ray onset radius and the radius for unit radial optical 
depth, i.e. Ri = max[_Ro, -Ri]. For a wind with mass loss 
rate M and flow speed Voo, this radius for transition from 
optically thick to thin X-ray emission is given by 



Ri 



T,R, = 



KhfM 



25Rp 



(23) 



where r* is a characteristic wind optical depth to the stel- 
lar surface radius i?*, and the latter equalities use the scal- 
ing ((S]) for the bound-free X-ray opacity Kbf, with M-e = 
M/(10-''Mq yr-i) and l/iooo = V^o/lOOOkm s'^ 

In direct analogy, we can similarly deflne a characteris- 
tic adiabatic radius for transition from radiative to adiabatic 
cooling of the associated wind shocks, 



Rm 



47rVcx 



625i?p 



M- 



■0 ■ 



(24) 



where Tkev ~ kTs/keV. For X-rays emitted with photon en- 
ergy comparable to the shock energy, -Bkov ~ Tkcv, we flnd 
Rsl/Ri « Kca/nhi ~ 25. As detailed below, this implies there 
can be extensive wind regions in which the density is high 
enough for shocks to be radiative, but too low for much wind 
absorption of their emitted X-rays. In particular, indepen- 
dent of the mass loss rate, any shocks formed near the radius 
Ri, where X-rays have near unit optical depth, should al- 
ways be strongly radiative. 

Figure [2] plots estimates of R^^/Ro vs. O-star spectral 
subtype, for shocks with the temperature Tkev = 0.5, as ex- 
pected for wind-shocks generate d by t h e line-deshadowing 
insta bility (|Runacres fc Owockil |2002| : iDessart fc Owockil 
l2003h . a nd estimated from the observed, relatively soft X-ray 
spectra (IWoidowski fc Schulj l2005l : IZhekov fc Pailal 120071 : 



iNaze et al.ll2012h . The large Ra. values conflrm that wind- 
shocks of this energy should be radiative. 

This plot can also be readily used to estimate expected 
unit optical depth radii. In particular, for an emission line 
at -Ekcv = 1, the corresponding Ri would be a factor 100 
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Figure 2. Ratio Rn/Ro of adiabatic radius to shock onset radius, 
assuming Ro = 1.5-R*, and a post-shock temperature Tj^ov = 0.5. 
To co nstruct the plot , stella r parameters were taken from Tables 
1-3 of iMartins et al.l ||2005| ') {T^s, Lbol, R*, M,) for luminosity 
classes V, III, and 1 (marked respective by circles, squares, and 
triangles). The terminal spe eds are computed as Voo = S.GVosc, 
and mass loss rates are from IVink et al. 1 1 I2OOOI) . 



smaller than the _Ra plotted in figure (2] Since this means 
most 0-stars would have Ri/Ro < 1, we see that absorption 
effects should be weak to moderate throughout the O-star 
domain, even while the shocks are generally radiative. See ^ 
and figure U for further discussion and illustration of these 
scalings. 



3.2 Bridging law with thin-shell mixing 

Previous analyses using the density-squared emissivity p6p 
(e.g., OC99, OCOl) have directly parameterized the X-ray 
filling factor /v as following some specified radial function, 
e.g., a power-law. But instead let us now parameterize the 
shock number distribution in (|17|) by an analogous power- 
law. 



Tisir) ; 



— V ; r>R. 
r 



(25) 



with Ws = for r < Rp. Both instability simu- 
lations (lOwocki fc Pu5 I1999I : iRunacres fc Owockil I2OO2I : 
iDessart fc Owockil booj) and X-ray profile fitting and He- 
like f/i rati os (|Leutenegger et all I2OO6I : ICohen et all |2006| . 
l20ld . [ 20 111 ) suggest an initial onset for shock formation 
around Ro ~ 1.5R,. 

Applying (|25p in (|17|l , the analysis in i)2.4l provides a 
more physical model of shock X-ray emission that accounts 
explicitly for the effects of radiative cooling. However, it still 
does not account for any thin-shell mixing. The inherent 
thinness {£ <C r) of radiative shock cooling zon es makes 
thern subject to various thin-shell instabilities (jVishniad 
1 19941 ). which in numerical simulatio ns lead to highly com - 
plex, turbulent shock str ucture (e.a.. lWalder fc Folini|[l99^ . 
iParkin fc PittardI (|2010l ) discuss how the inherently limited 
spatial resolution of radiative shock simulations leads to a 
"numerical conduction" that transports heat from liigli to 
low temperature gas, resulting in a severe, but difficult-to- 
quantify reduction in the X-ray emission. 



While the specific mechanisms within hydrodynamical 
simulations may indeed depend on such numerical artifacts, 
the perspective advocated here is that such an overall reduc- 
tion in X-ray emission is likely a natural consequence of the 
turbulent cascade induced by the thin-shell instability; this 
should lead to substantial physical mixing between cool and 
hot material, with the softer and more efficient radiation of 
the cooler gas effectively reducing the emission in the X-ray 
bandpass. 

Pending further simulation studies to quantify such 
mixing and X-ray reduction, we make here the plausible 
ansatz that the reduction should, for shocks in the radiative 
limit l/r <^ 1, scale as some power m of the cooling length 
ratio, {l/ry^'. To ensure that the mixing becomes ineffective 
in the adiabatic limit - for which the cooling layer is too 
extended to be subject to thin-shell instability - we recast 
the bridging law pO|l for the shock volume filling factor /vs 
in the generalized form, 

~ (l + r/^)i+™ " (1 + Kcpr)i+™ ' ^^^^ 

with the level of mixing now controlled by a positive value of 
the mixing exponent m. To simplify notation in the analysis 
to follow, we have dropped here the subscripts ( "s" or "w" ) 
for the quantities (e.g., £, r, Kc, p) on the right-hand-side. 
Recalling that the cooling coefficient takes the scalings 
given in equations H8a|) - (|8c)) . we can use the adiabatic radius 
7?a from (|24|) to characterize the asymptotic regimes for (|26|) . 

If i?a < ^^o, then the shocks are adiabatic through- 
out the wind, and we recover the standard density-squared 
emissivity (|16|l with a specified volume filling factor set by 
/v = ris. 

If -Ra > -Ro, this adiabatic scaling still applies in the 
outer wind, r > Ra,; but in the inner regions Ro < r < Ra_, 
the cooling term in the denominator of (|26p dominates, giv- 
ing the emissivity (|16|) now a reduced dependence on density. 



'7x 



i?o < r < i?a 



(27) 



' (Kcr)i+™' 

Without mixing (m = 0), the density dependence is thus 
linear, but with mixing (m > 0), it becomes sub-linear. 

As noted above, absorption is a modest effect in even 
dense O-star winds, with r« at most of order unity, implying 
then that Ri < Ro (jCohen et al.ll2O10h . But the stronger co- 
efficient {nc/fibi = Ra/Ri ~ 25) for radiative cooling means 
that such winds typically have i?a > Ro , implying that most 
shocks remain radiative well above the wind acceleration re- 
gion where they are generated. 

3.3 Lx scalings with M /Voa 

Let us now turn to the scalings for the overall X-ray lumi- 
nosity. Applying the emissivity (|16|l and filling factor p6|l 
to the exospheric model (|22p with the power-law shock dis- 
tribution (|25p . we find 
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(1-f Kcpr)i+'^ 



■ r^dr 



dr 



rP{rwY-^{rw + RaY+^ ' 

CsUsoRo, and the latter equality uses mass con- 
servation to cast the integral in terms of the scaled wind 



where Cn 
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velocity, w{r) = v{r)/Voo- This is generally taken to have 
a 'beta-law' form w{r) = (1 — R,/r)^ , with the canonical 
case /3 = 1 giving rw = r — R,, and the constant-speed case 
/3 = giving rw = r. For these special cases, ij5] gives some 
results for full integrations of (|28|) . 

But even for a generic velocity law w(r), we can read- 
ily glean the essential scalings of the Lx from (|28p in key 
asymptotic limits of adiabatic vs. radiative shocks. 



3.3.1 Adiabatic shocks in optically thin wind 

First for low-density winds with optically thin emission from 
adiabatic shocks, Ri < Ra, < Ro, we can effectively drop the 
Ra, term in the denominator, and set the lower bound of the 
integral to the fixed onset radius, _Ri = Ro, yielding 



■ M' 


2 


r°° dr 









Rl < i?a < Ro ■ 



(29) 



Since the resulting integral then is just a fixed constant that 
is independent of M, the X-ray luminosity recovers the stan- 
dard adiabatic scaling Lx ~ (M/Voo)^. 



3.3.3 Summary of asymptotic scalings 

To summarize, power-law shock-number models with thin- 
shell mixing have the asymptotic scalings, 

-I 2 



M 



M 



M 



; Rl < Ra < Ro ; adiabatic, thin (31a) 

1 — m 

; Rl < Ro < Ra, ; radiative, thin (31b) 

i-p 

; Ro < Rl < Ra ; radiative, thick (31c) 



where the progression represents a trend of increasing 
M/Voo. The first applies for weak winds, for example from 
early B main-sequence stars, which typically s how weak X- 
ray e mission, with Lx well below 10~'^Z/boi (jCohen et al.l 
The middle scaling for intermediate-density winds is 
the most relevant for the bulk of O-type stars found to follow 
the Lx ~ Lboi relation. The last applies only to the strongest 
winds, e.g. very early O supergiants like HD93129A, for 
which analysis of X-ray li ne profiles show mo derate absorp- 
tion eflfectCl, with r, > 1 (|Cohen et al.ll201ll '). 



3.3.2 Radiative shocks m optically thin or thick wind 

For high-density winds with radiative shocks and so Ro < 
Ra, the Ra now dominates its term, and so can be pulled 
outside the integral. Rescaling the remaining integrand in 
terms of the initial radius _Ri, we find 



Lx — Co 



M 



Vex 



Ra. 



dr 



rP{r'w)^' 



Ro < Ra, (30) 



where again the integral is now essentially independent of 
M/Voo. 

For the intermediate-density case in which the radiative 
shock emission is optically thin, Ri < Ro < Ra,, the inte- 
gration lower limit is fixed at the onset radius, Ri — Ro, 
which is independent of M/Voo. But since Ra ~ M/Voo, the 
overall scaling is Lx ~ (M/Voo)^""". 

For the case of very dense, optically thick winds with 
radiative shocks, Ro < Ri < Ra, the lower boundary at 
-Ri — Rl gives the residual integral an additional depen- 
dence on _R™~^; since Ri too scales with mass loss rate, the 
dependence on the mixing index m cancels. At the radius 
Rl the cooling length ratio £/r is always the same, implying 
that in optically thick winds the observed radiative shock 
emission is likewise fixed for any mass loss rate. The lumi- 
nosity scaling thus becomes independent of mixing, scaling 
just with shock number index as Lx ~ (M/Voo)^"''. 

This is in fact the same scaling found in OC99 for opti- 
cally thick winds, but with adiabatic shocks. Indeed, OC99 



argued that assuming a filling factor /v 



could give 



a sub- linear dependence on mass loss rate, Lx ~ M"'®, and 
so possibly reproduce the Lx ~ Lboi relation. Subsequent 
analysis of X-ray line profiles observed from Chandra and 
XMM-Newton have shown, however, that optical depth ef- 
fects are quite mod erate even for O-st ars like Puppis with 
quite dense winds (|Cohen et al.|[201ol '). The bulk of O-star 
winds are simply too low-density for this optical thickness 
scaling to apply, and so this cannot be the explanation for 
the Lx ~ Lboi relation. 



3.4 Link between M and Lboi scaling 

As noted in the introduction, and summarized in figure [T] 
straightforward application of CAK wind theory implies a 
direct dependence of mass loss rate on luminosity that scales 
as M ~ L\I^ , where a « 0.6 is the CAK power index; for 
the bulk of O-type stars with intermediate-density winds 
and thus Lx ~ AL^~'", reproducing the observed Lx ~ Lboi 
relation thus simply requires a mixing exponent m ~ 1 — a ~ 
0.4. 

More generally, let us now consider how this require- 
ment is affected if one accounts also for a secondary depen- 
dence on stellar mass M, which in turn can give a further 
indirect dependence on Lboi. Specifically, within CAK wind 
theory, M ~ M^"^/" and 14o ~ M^^^, and so if we in turn 
assume from stellar structure a mass- luminosity dependence 
M ~ L;;„i, where s « 1/3 (e.g.. [Maeder,,200g. p. 360), we 
find 



M\ 2-s(2~a) log(Lx) 
log k log Lboi ~ -r^ — - . (32) 



where the latter relation makes use of the scaling H31bp 

Reproducing the empirical Lx ~ Lboi relation thus now 
requires 

_ 2(l-a)-s(2-a) 



2 - s(2 - a) 



(33) 



Accounting for a stellar structure scaling s ~ 1/3 with a 
CAK index a « 0.6 thus now requires a mixing exponent 

^ Indeed, this star could be viewed as a transitional object to the 
WNH type Wolf-Rayet stars, for which absorption effects should 
strongly attenuate any X-rays from instabilities in the wind ac- 
celeration region; see 

* Of course, the terminal speed also depends on stellar radius as 
Voo ~ 1/ V Rt, , but the diverse luminosity classes of X-ray emit- 
ting O-stars makes it difficult to identify any systematic depen- 
dence on Lboi that might influence the overall Lx — Lboi relation. 
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Figure 3. X-ray line profiles (normalized to unit maximum) 
vs. Doppler-shifted wavelength (normalized to shift for terminal 
velocity Voo), for optically thin (r* <SC 1; left) and marginally 
optically thick (t» = 1; right) lines. The black curves show the 
adiabatic, constant-filling-factor models that give good general 
fits to observed X-ray emission lines. The other curves show ra- 
diative models for no mixing and for m = 0.4, with shock-number 
indices p as labeled. Radiative models with constant shock num- 
ber (p = 0) give profiles that are too broad (red curves). But 
models with steeper number exponents (p = 1.5 for m = in 
green; and p = 2 for m = 0.4 in blue) fit well the observationally 
favored black curves. If the O-star Lx — ibol relation is to be re- 
produced with moderate thin-shell mixing (m 0.4) of radiative 
shocks, then matching observed X-ray emission lines requires the 
shock number to have a moderately steep decline (p 2) above 
the initial onset radius R^. 



of radiative cool i ng aff ects X-ray line profiles. Following 
lOwocki fc CohenI (|200ll . hereafter OCOl), the directional 
Doppler-shift of the X-ray line emission within the expand- 
ing wind is modeled through a line emissivity r]x{r,^) at 
an observer's wavelength A along direction cosine p from 
a radius r. The resulting X-ray luminosity spectrum Lx is 
computed from integrals of the emission over direction and 
radius, attenuated by bound-free absorption within the wind 
(cf. OCOl equation 1), 

/I /■ oo 

dii / drr2j7A(M,r)e""'^'''^ (34) 
-1 J 

The absorption optical depth r(^, r) is evaluated by convert- 
ing to ray coordinates and then integrating for each ray with 
a fixed impact parameter from the local position to the ob- 
server. For the standard P — 1 velocity law, the integrals are 
analytic, with overall scaling in proportion to r, = Ri/R,. 

In principle, this integrated optical depth can be af- 
fected by the "porosity" associate d with optically thic k 
clumps or anisotropic "pancakes" (|Feldmeier et al.l |2003| ) , 
with potential consequences for interpreting the asymme- 
try of X-ray line profi les in terms of the wind mass loss rate 
(|Oskinova et al.ll2006l ) . But for the optically thin (r, < 1) or 
marginally optically thick (r« ~ 1) lines c onsidered here, in- 
dividual clumps should be optically thin (jOwocki fc CohenI 
l2006l : ISundavist et~al]|2012l ) , and so such effects are not im- 
portant for the discussion below, which focuses on the overall 
profile width. 



m « 0.22, somewhat smaller than the m ~ 0.4 required 
if one assumes no systematic mass- luminosity scaling (i.e., 
s = 0). 



4 EFFECT OF SHOCK COOLING ON X-RAY 
LINE PROFILES 

4.1 Basic formalism 

Observations by XMM-Newton and Chandra of spectrally re- 
solved, wind-broadened X-ray emission lines from luminous 
OB stars provide a key diagnostic of the spatial distribu- 
tion of X-ray emission and absorption in their expandin; 
stellar winds ([ignace fc Gavlevl[2003 : lOskinova etahlbOOL. 
lOwocki fc CohenI 120061 " In particular, the relatively weak 



blue-red asymmetry of the observed emission lines indicates 
modest wind optical depths, r» ~ 1, while the overall width 
constrains the spatial location of the emission within the 
expanding velocity law. For the usual density-squared emis- 
sion model with a prescribed (power-law) spatial variation in 
volume filling factor /v , fits to observed X-ray lines are typ- 
ically consistent with a standard /3 ~ 1 velocity law and an 
/v that is spatially nearly constant, corresponding to p « 
within th e adiabatic scaling implicit in the emiss ion measure 
analysis (|Kramer et al.ll2003l : [Cohen et al.ll2006l ). In the dis- 
cussion below, we refer to these profiles - plotted in black 
in figure 13]- as "observationally favored". 

Within the perspective discussed here that shocks 
within most O-star winds are likely to be radiative in- 
stead of adiabatic, let us now examine how inclusion 



4.2 Scaling analysis 

As noted, applications of this OCOl formalism assuming a 
density-squared emission show that a spatially constant X- 
ray volume filling factor /v gives generally quite good fits to 
observed X-ray lines. To examine the effect of radiative cool- 
ing, let us now apply the more general bridging-law scalings 
of p6p and (|26|l . assuming the simple power-law form for 
the shock number p5|l . The spatial integration thus takes 
the same form as the exospheric result (|28p . except that ab- 
sorption is now treated explicitly by the exponential optical 
depth term in the integrand, with the radial lower bound 
fixed at the X-ray onset radius, Ri = Ro- 

We can again infer basic scaling results by inspection of 
this integrand in the limits of adiabatic vs. radiative shocks. 
For the adiabatic case, this follows the scaling in (|29|l . with 
radial dependence To match observed profiles, 

such adiabatic emission models require constant /v = ris, 
with the zero power-law exponent p — implying a l/{wr)^ 
variation of the integrand. 

By contrast, the radiative limit follows the scaling in 
(|30|l . with direct radius dependence l/r^^^""", and velocity 
dependence l/w^~"'' that is weaker than the of the 

adiabatic model. Focusing first just on the former, we see 
that reproducing the 1/r^ integrand needed to fit observed 
profiles would now require p = 1 -I- m. In practice, to com- 
pensate for the weaker inverse-speed dependence, fitting the 
observed profile width requires a somewhat steeper shock- 
number decline, as we now quantify. 
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4.3 X-ray line profiles from radiative shocks 

For a /3 = 1 velocity law with Ro — 1.5R* and various 
specified values of the exponents p and m, figure [3] plots 
normalized X-ray line profiles for optically thin (r, ^ 1; left) 
and marginally optically thick (r* = 1; right) lines. In both 
cases, the black curve represents the adiabatic, constant- 
filling-factor (p = 0) model that gives generally good fits to 
observed profiles. 

The other curves show results for radiative shocks. 
Without mixing, the red curve for radiative shocks with con- 
stant shock number {p = 0) is far too broad; fitting the fa- 
vored black profile now requires up — 1.5 (blue curve), which 
is even steeper than the predicted p = m -f 1 = 1 needed to 
compensate for the weaker direct radial scaling. With mixing 
exponent of m = 0.4, we find a p = 1.5 shock-number model 
gives profiles (not shown) that are still somewhat too broad. 
But with a somewhat steeper number exponent p = 2, the 
blue curve again nearly reproduces the black curviQ. 

An overall conclusion is thus that, for radiative shock 
models with mixing at a level needed to reproduce the Lx ~ 
Lboi relation, matching observed X-ray emission requires a 
steep radial decline (p ~ 2) in shock number above the onset 
radius. 



4.4 Decline of Lx in optically thick winds 

Such a steep radial drop off in shock number has important 
implications for the scaling of X-ray luminosity for the high- 
est density stars that become optically thick to bound-free 
absorption, i.e. with r, > 1 and thus _Ri > Ro- Namely, 
from the scaling given by H31cp . we see that taking p ~ 2 
implies that the X-ray luminosity for such optically thick 
winds should now decline inversely with mass loss rate, 
Lx ~ M^~^ ~ 1/ M. If the X-ray emission is concentrated 
near an onset radius within the wind acceleration zone, the 
bound-free absorption by the overlying, optically thick wind 
significantly attenuates the net X-rays seen by an external 
observer. For early 0-type supergiants with dense winds, 
this can lead to a reduced X-ray luminosity, but because the 
overall decline of bound-free opacity with X-ray energy, the 
observed spectrum can be hardened. The recent analysis of 
X-rays from the 02If star HD9 3129A provides a potential 
example approaching this limit (|Cohen et al.ll201ll '). 



5 SCALING RESULTS FOR FULL 
EVALUATION OF Lx INTEGRAL 

5.1 Constant-speed wind with /? = 

As a supplement to the asymptotic Lx scalings given in i]3.3l 
let us finally consider full solutions for the general integral 
((28]). For a constant speed model with (3 = Q and thus wir) = 
1, general analytic integration is possible in terms of the 
incomplete beta function. But the typical properties can be 



Also not shown here are profiles computed for p = 2 and the 
alternative mixing exponent value m = 0.22, which we find also 
give close agreement with the black curves. 



more simply gleaned by examining the special case p = 1, for 
which the integral in (|28|) takes the simpler analytic form, 



Lx = 



1 + mRJRi 



1 



(35) 



167r2m(l-m) [{1 + R^/Ri) 

wherein the square-bracket factor sets the scalings with 
M/Voo, with the preceding terms just fixing the overall nor- 
malization. For low-density, optically thin winds, the ini- 
tial radius is fixed to the onset radius, so that R^/Ri ~ 
Rb./Ro ~ M/Voo. For Ra/Ro <^ 1, expansion of the square- 
bracket term recovers the adiabatic scaling Lx ~ (M/Voo)^ 
of piafl . while for R^/Ro ^ 1, it becomes proportional 
to Ri^"^ and so recovers the radiative, optically thin scal- 
ing (|31b|) . For high-density, optically thick winds, the ratio 
Ra/Ri = Ra/Ri becomes independent of mass loss rate, and 
so the square-bracket factor and thus Lx approach a con- 
stant value, in accord with H31cP for this case with p = 1. 



5.2 Standard wind with /3 = 1 

For a standard /3 = 1 velocity law with rw = r — Rt,, ana- 
lytic integration of (|28|l can be cast in terms of the Appell 
hypergeometric function; but in practice it is more straight- 
forward just to carry out the integration numerically. 

Figure m plots the resulting X-ray luminosity Lx (scaled 
by the dimensional factor outside the square brackets in 
(|35|l ) vs. adiabatic radius 7?a (scaled by the shock-onset ra- 
dius -Ro), as computed from numerical integration of (|28|l 
for the selected, labelled values of the mixing exponent m 
and shock number exponent p. As seen from (|24|l . plotting 
vs. Ra. can be viewed as a proxy for plotting vs. the wind 
density parameter M/Voc- 

The left vertical line at R^/ Ro ~ 1 represents the transi- 
tion from low-density adiabatic shocks, for which Lx ~ M'^, 
to intermediate-density radiative shocks, for which Lx ~ 

The right vertical line represents the transition from 
optically thin to thick winds, implemented here through 
a sudden change in the integration lower bound, Ri — 
max(_Ri, i?o), with the shock onset radius fixed at Ro — 
1.5J?«, and Ri the radius for unit optical depth. Like _Ra, 
Til scales with A'l/Voo, with values that are a small, fixed 
fraction Ri/Rn = fia of the adiabatic radius. Since ob- 
served X-rays are typically a modest factor higher energy 
than the characteristic wind shock (e.g., iSkev « 1 « 2rkcv), 
the curves plotted in figure 3] assume, following (|23|l and 
dull, /la = i?i/-Ra = Kbf/Kc ~ 1/100. For Very large 
-Ra > Ro/fia ~ 100-Ro, we thus have Ri > Ro, leading 
to a declining Lx, as predicted by the optically thick wind 
scaling Lx ~ Ri~'' from (|31c|l . 

But for moderately dense winds, with Ro < Ra < 
Ro/fia (between the vertical lines in the figure), the in- 
creasing Lx approaches the power-law variation Ri~"^ pre- 
dicted by (|31b|) . In particular, the black curve with p = 2 
and m = 0.4 represents the preferred model with sub-linear 
scaling in J?a, and thus in M/Voo, implying a nearly linear 
scaling of Lx with Lboi- 

Specifically, for typical values for stellar radius (i?* ~ 
10 — 2O-R0) and wind terminal speed {Voo ~ 2000 km s^^), 
this intermediate-density regime with Lx ~ Lboi applies to 
wind mass loss rates that range from below lO~^Af0 yr~^ to 
a few times lO~^M0yr~^; this essentially encompasses the 
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Figure 4. Normalized X-ray luminosity Lx vs. adiabatic radius i?a scaled by shock onset radius Ro = 1.5-R*, for winds with a standard 
/3 = 1 velocity law, and selected, labelled values for the parameters p and m . The dashed, dot-dashed, and dotted lines show the expected 
scalings for respectively adiabatic, radiative, and thin-shcU-mixed shocks. As discussed in the text, the variation in represents a proxy 
for the wind density parameter M/Voo, ranging from low density, adiabatic shocks on the far left, to high-density, optically thick winds 
with radiative shocks on the far right; the optically thick turnover at Rs^/Ro > 100 applies for observed X-ray energies that are about 
a factor two higher than the shock energy. The intermediate density case with optically thin, radiative shocks follows the Lx ~ M^^*" 
scaling that reproduces the empirical Lx ~ Lboi relation if m, ft! 0.4. The upper axis uses the analysis from figure [2] to mark the 
corresponding spectral class for main sequence stars; for higher luminosity giants and supergiants, the associated spectral class sequence 
would shift to the right. 



entire O-star spectral range for which the Lx ~ 10 ^Lboi 
relation is found to hold. 



6 CONCLUDING SUMMARY 

The central result of this paper is that, in the common 
case of moderately dense O-star winds with radiative shocks 
(i?a > ^o), thin-shell mixing can lead to this sub-linear 
scaling of the X-ray luminosity with the mass-loss rate, 
Lx ~ (M/V^tx,)^"™. Depending on the secondary scalings 
of wind density with bolometric luminosity, one finds that 
m ~ 0.2 — 0.4 can give roughly the linear Lx-I/boi law that is 
empirically observed for O-star X-rays. Further simulation 
work will be needed to see if such mixing exponent values 
are appropriate, and indeed to test the validity of the basic 
mixing exponent ansatz. 

But in the course of exploring this idea of thin shell 
mixing, the analysis here has lead to several interesting sec- 
ondary results with validity and implications that are largely 
independent of mixing or any specific model for it. A sum- 
mary list includes: 

(i) In contrast to previous analyses that invoked a 
density-squared emission measure description for shock pro- 
duction of X-rays, we derive here a more general bridging 



law showing how the density-squared scaling of adiabatic 
shocks transitions to a single-density scaling for radiative 
shocks. 

(ii) For radiative shocks, the X-ray volume filling factor 
is not fixed (as is commonly assumed), but is reduced by the 
narrow extent of the shock layer, /v ~ £/r. 

(iii) For nearly all O-stars, the large radiative-adiabatic 
transition radius, Ra. ^ Ro ~ 1.5/?,, implies that 
instability-generated shocks in the wind acceleration region 
should follow the radiative scaling, giving the X-ray lumi- 
nosity a linear scaling with mass loss rate, Lx ~ M/Vaa- 

(iv) For low-density winds of lower luminosity (early B) 
stars, shocks should indeed become adiabatic, implying then 
a steep (Lx ~ M^) decline of X-ray luminosity, as is in 
fact generally found for single, non-magnetic, early B-type 
stars, for which the inferred X-r ay emission measu r e ofte n 
approaches that of the full wind jCohen et al.lll997l . 120081 ') 

(v) Matching observed X-ray emission lines with such 
models of radiative shocks with or without thin-shell mixing 
requires the shock number to have a moderately steep de- 
cline above the X-ray onset radius, with power-law exponent 

1.5-2. 

(vi) This in turn implies that the scaling of X-ray lumi- 
nosity for dense, optically thick winds should saturate and 
even decline with increasing mass loss rate. 
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This last result on X-ray absorption is perhaps not too 
relevant for most of the 0-stars following the Lx — Lboi re- 
lation, for which optical depth effects are weak to marginal. 
But it can become important for the dense, moderately op- 
tically thick winds of extreme, eaxly O-stars like HD93129A, 
which can be viewed as a transitional object to the WNH 
type Wolf-Rayet stars. More generally, the high-density of 
Wolf-Rayet winds imply that absorption should strongly at- 
tenuate X-rays from any instability-generated shocks in their 
wind acceleration region. As such, the observed hard X-rays 
seen from Wolf-Rayet stars like WR6 (EZ CMa) seem un- 
likely to be explained by this standard model of LDI shocks 
(|Qskinova et al.ll2012l '). This also has potential implications 
for interpreting observed X-rays from very massive stars that 
have Lx ~ 10~^I/boi despite ha ving very high wind opti- 
cal depths jCrowther et aLlbOld l. and whether these might 
instead originate from wind- wind collisions of close, unde- 
tected binary companions. 

Indeed, the mixing ansatz in this paper could also be ap- 
plied to model X-ray emission from colliding wind binaries, 
and their Lx scaling with orbital separation. Wide binaries 
with adiabatic shocks should still follow the usual inverse 
distance scaling, as directly confirmed by observations of 
mult i-year-period ell iptical systems like WR140 and 77 Cari- 
nae (|Corcoranll201ll ). But in close, s hort (day to week) pe - 
riod binaries with radiative shocks (|Antokhin et al.ll2004h . 
mixing could reduce and li mit the effective X-ray emis- 
sion from the wind collision (|Parkin fc Pittardll2010l '). and 
thus help clarify why such systems often ha rdly exceed the 
Lx ~ 10~^Lboi scaling found for single stars (|Oskinovall2005l : 
ICorcoran|[201ll : lGagnell201ll ). 

Finally, in addition to exploring such effects in colliding 
wind binaries, a top priority for future work should be to 
carry out detailed simulations of the general effect of thin- 
shell mixing on X-ray emission, and specifically to examine 
the validity of this mixing-exponent ansatz for modeling the 
resulting scalings for X-ray luminosity. 
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